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ABSTRACT 

Detailed information on solar radiation characteristics on Mars are necessary for effective design 
of future planned solar energy systems operating on the surface of Mars. In this paper we present a 
procedure and solar radiation related data from which the daily variation of the global, direct beam and 
diffuse insolation on Mars are calculated. Given the optical depth of the Mars atmosphere, the global 
radiation is calculated from the normalized net flux function based on multiple wavelength and 
multiple scattering of the solar radiation. The direct beam was derived from the optical depth using 
Beer’s law, and the diffuse component was obtained from the difference of the global and the direct 
beam radiation. The optical depths of the Mars atmosphere were derived from images taken of the 
Sun with a special diode on the cameras used on the two Viking Landers. 

1. INTRODUCTION 

Missions to the Martian surface will require electric power. Of the several possibilities, 
photovoltaic power systems can offer many advantages, including high power to weight, modularity, 
scalability, and a long history of successful application in space. To design a photovoltaic system, 
detailed information on solar radiation data on the Martian surface is necessary to allow more accurate 
estimates of photovoltaic power system size and mass in system analysis and trade-off studies of 
relevant technology options. 
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Solar irradiance on the Mars surface varies with season and with the amount of atmospheric 
dust. A major concern is the dust storms, which occur on both local and global scales, and their effect 
on solar array output. Global storms may persist for long periods of time, such that the requirement 
for energy storage quickly becomes much too large to be practical. However, as we have shown there 
is still an appreciable large diffuse component of sunlight available even at high optical depths of the 
Mars atmosphere, so that solar array operation is still possible. Mars may be considered "clear” most 
of the time for most places, but when the global storms occur, their effects are large. 

The variation of the solar radiation on the Martian surface is governed by three factors: ( 1 ) the 
variation in Mars-Sun distance, (2) variation in solar zenith angle due to Martian season and time of 
day, and (3) opacity of the Martian atmosphere. Opacity measurements were derived from images 
taken of the Sun with a special diode on the cameras used on the two Viking Landers. Viking 
Lander 1 (VL1) is located at 22.3° N latitude and 47.9° W longitude (in Chryse) and Viking Lander 2 
(VL2) is located at 47.7° N latitude and 225.7° W longitude (in Utopia). Mars maps can be found in 
[1]. The camera observations have produced 1044 measurements of optical depth spread over 1-1/3 
Martian years. 

The global solar radiation is composed of two components. The direct (or "beam”) component 
is the light which comes directly from the Sun, attenuated by atmospheric dust. The diffuse 
component is the light reaching the surface from atmospheric scattering. Given the optical opacity, 
global radiation is calculated from the normalized net solar flux function based on multiple wavelength 
and multiple scattering of the solar radiation. The direct beam component was derived from the 
optical depth of the atmosphere using Beer’s law, and the diffuse component was obtained from the 
difference of the global and direct radiation. 

The Martian atmosphere consists mainly of suspended dust particles, the amount of which vary 
daily, seasonally, and annually, depending on local and global storm intensities and their duration. 

The intensity of the dust storms are defined in terms of the optical depth of the Martian atmosphere. 
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From a photovoltaic system design point of view, the intensity, frequency, and duration of these 
storms may be viewed as "partially cloudy” and "cloudy” days for which additional energy storage 
must be taken into account. 

•>.» This paper is an update to reference [2] and includes a refinement of the solar radiation model; 
Mars surface albedo; optical depths at different latitudes; and solar radiation data for the entire planet. 
The results are presented in a series of figures and tables. The solar radiation data and the procedure 
presented in this paper may be used for the calculation of any desired solar radiation quantity on Mars 
in an engineering design. 

2. OPTICAL DEPTH AT VIKING LANDERS 

The most direct and probably most reliable estimates of opacity are those deri ved from the 
Viking Lander imaging of the Sun, Figures 1 and 2 show the seasonal variation of the normal 
incidence of the optical depth at the Viking Lander locations VL1 and VL2, respectively. The season 
is indicated by the value of L s areocentric longitude of the Sun, measured in the orbital plane of the 
planet from its vernal equinox, L s = 0°. Figures 1 and 2 were derived from [3] and were discretized 
for each 5° (see measured and discretized). Opacities are minimum during the northern spring 
(L s = 0® to 90°) and summer (L s * 90° to 180°), and maximum during southern spring (L s = 180® 
to 270°) and summer (L s = 270° to 360°), the seasons during which most local and major dust storms 
occur. When dust storms are not present, the optical depth is typically about 0.5. 

Two global dust storms occurred during the periods of each observation, as indicated by the 
high values of the optical depth (they are lower bound values). Numerous local dust storms were 
observed as well. Figures 3 and 4 show the variation of the optical depth of Figs. 1 and 2 in greater 
detail. 

The cameras on the Viking Lander measured the optical depth in the morning and in the 
afternoon. The morning optical depths are higher, on the average, than the afternoon optical depths. 
These differences indicate a diurnal cycle of water condensation and evaporation. The afternoon 
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optical values are more representative of that solely due to the suspended dust particles. These values 
are shown in Figs. 1 and 2 and are assumed to remain constant throughout the day for the solar 
radiation calculation. Other methods [4,5] for the estimation of the optical depth of the Mars 
atmosphere deal with die solar reflected or thermally emitted radiation. The optical depth derived 
from the Viking Lander imaging of the Sun, which is used in the paper, is a direct method and is 
independent of the optical properties and temperatures of the ground and airborne dust. 

3. OPTICAL DEPTH AT DIFFERENT LATITUDES 

The optical depths at different latitudes other than the ones measured at the Viking Landers 
VL1 and VL2 are introduced in this article and are based on the assumption that the optical depth 
varies spatially and are expected to be greatest in the source region where dust is raised into the 
atmosphere. References [4,5] support this assumption. It is therefore assumed that the optical depths 
have a Gaussian distribution centered at the times of the global storms at VL1 and VL2. The 
relationship of the optical depth with latitudes may be expressed by a low order polynomial ratio. The 
optical depth function coefficients correspond to the measured opacities at VL1 and VL2. This 
function: t(<J), L s ), was developed for the case of two global dust storms per Martian year; a model 
that corresponds to the Viking Lander observation. 

The location of the first global storm is well characterized, however there is some uncertainty 
about the location where the second 1977 storm started. We have therefore developed two models for 
global dust storm opacity, corresponding to latitude -10° and -30° locations for the second storm. The 
primary source areas for dust during the growth phases of global dust storms on Mars have historically 
been the Hellas and the Solis Planum regions. For the first model we thus assumed that the two 
global dust storms have originated both at latitude <j> = -30° and took place and at L sl = 215° (first 
storm) and at L s2 = 295° (second storm). The model for the optical depth function, is assumed to 
have the form: 
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i.e., the optical depth takes the greater value between 0.5 and the value in the square parenthesis. The 
numerical constants were determined by several conditions based on the Viking Lander measurements 
and least square fitting. Figures 5 and 6 describe the variation of the optical depth with latitude <J> 
and areocentric longitude L s based on Eq. (1). 

Reference [6] states that the second 1977 storm was centered in the Valles Marineris area 
( 10° S, 50° to 80° W). A second optical depth function was therefore developed with the first global 
storm at latitude <J>j = -30° and the second global storm at latitude <j> 2 “ -10°. The times of the 
storms were taken the same as before at L s = 215° and L s = 295°. In this case, the optical depth 
function was assumed to have the following form: 


T(<j>X,) = maxi 0.5, 19,500 



-(L -215)7730 


4,000 + (4> + 48. ir 


♦ JL| 

+ 12,700 410 i 

2,465 + (<J> + 13.1 ) 2 

i.e., the optical depth takes the greater value between 0.5 and the value in the square parenthesis. The 
numerical constants were determined similar to the preceding model. Figures 7 and 8 describe the 
variation of the optical depth with latitude <|> and areocentric longitude L s based on Eq. (2). 

4. DUST STORMS 


Observation from Earth, from orbit around Mars, and from the surface of the planet itself have 
shown that airborne dust is a highly variable and thermodynamically significant component of the 
Martian atmosphere. Martian dust plays a key role in determining the current climate of Mars, and is 
suspected of having a major influence on the evolution of the surface and the history of climatic 
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conditions on the planet. Airborne dust may be transported over short or long distance depending 
upon regional meteorological conditions. Dust is raised from the surface into the atmosphere, and 
when certain conditions are met, it may be translated into local or global dust storms. The spatial 
distribution of dust storms and of changes in surface albedo suggest that dust can be removed from the 
surface almost anywhere on the planet. There may be long periods when little dust moves at all at 
some locations, while there appear to be other regions where dust is raised more frequently. Many of 
the parameters controlling the location, severity, and frequency of dust-raising activity are not well 
understood. The classical view is that dust storm activity is correlated with a general intensification of 
winds around the time of maximum insolation, near southern summer solstice and perihelion. Dust 
storms are distinguished as minor or major storms, the former being localized clouds that remain 
locally and dissipate in a few days or so. Major dust storms are also called “global” or "great” 
storms and are classified as regional (do not encircle the planet) or planet-encircling (cover much of 
one or both hemispheres) storms. Local storms appear to play an important role in initiating the 
global storms. The thermodynamics of local storm generation and the mechanisms by which local 
storms become planet-encircling storms remain obscure. During the planet-encircling dust storms, dust 
is raised and spread over much of the planet for many Martian days (sols), and eventually falls to the 
surface. The direction of the dust transport of these storms is that dust is moved about within the 
southern hemisphere or from the south to the north over several weeks. The atmospheric dust 
decreases the amount of insolation reaching the surface, increases the infrared atmosphere heat flux, 
and absorb radiation emitted by the surface. A compilation record of all known clouds and dust 
activities on Mars reported in the literature over the past 1 15 years is reported in [7]. The 
characteristics of the global and local dust storms are discussed in the following subsections. 

4.1 Global Dust Storms (Regional and Planet-Encircling Storms) 

Global dust storms have been observed on Mars for many years, both from Earth-based 

. > " , v t: :■ . ... „ -'"V. " ■ ■ . “ ■■ ■ ■ . - ' ■< - ' 

telescope observations and from Mariner 9 and Viking spacecraft observations. These storms typically 
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develop during the southern summer season, presumably because wind speeds reach their peak at that 
time as the solar insolation is maximum at perihelion passage. The global dust storms appear to 
evolve from local storms in preferred regions. These regions - Claritas Fossae/Solis Planum, 
Hellespontus/Noachis, Sirtis Major Isidis (Fig. 9) - are characterized by their location in low and 
subtropical latitudes, by the presence of large east-west slopes, and by strong gradients in surficial 
albedo or thermal inertia. Global dust storms do not occur in northern spring and summer because of 
the greatly reduced insolation available at that time. Effects of major dust storms on the global 
opacity vary depending on the intensity and latitudinal extent of the storm. Initial buildup in opacity 
occurs over 2 weeks, decay begins immediately thereafter. The decay period of several months occurs 
between storms before a new storm is initiated. Although the major storms affected most of the 
planet, the latitudinal distribution is nonuniform. One or more dust storms of regional or larger size 
can occur in any given Mars year. There are many Mars years in which no planet-encircling dust 
storm has occurred; less frequently, there may be years in which not even regional storms occur. The 
duration of the global dust storms may vary from 35 to 70 days or more. The first global dust storm 
observed by Viking Landers (1977) spread from a latitude of 40 9 S to a latitude of 48° N in about 5 to 
6 days. The first storm moved northward at a speed of 10 m/sec and the second storm moved 
northward at a speed of 30 m/sec. Figure 1 also shows the seasonal date L s Of those Martian great 
dust storms (identified by the year in which they occurred) whose regional onset was obtained from 
Earth and/or by spacecraft or by the Viking Landers. When more than one storm has occurred in a 
single (terrestrial year), they are designated a or b (as in 1977a) in the order of their appearance. As is 
noticed, the storms were observed in southern spring and summer [8]. 

4.2 Local Dust Storms 

Three independent sources of information about local dust storms on Mars assess the spatial and 
temporal distribution of the clouds: (1) Viking Orbiter Infrared Thermal Mapper (IRTM); (2) Viking 
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Orbiter VIS (Visual Imaging Subsystem) data; and (3) the historical record of dust clouds occurrence 
produced from Earth-based telescopic observations [9]. 

Local dust storms, which may affect regions encompassing up to 10 km , have been observed 
on Mars during all seasons. However, they have been observed to occur most frequently in the 
approximate latitude belts 10° to 20° N and 20° to 40° S, with more clouds seen in the south than in 
the north and with more frequent sighting during southern spring and summer. Thus local storms are 
apparently most likely to occur during the same period as the global dust storms. The local storms are 
not uniformly distributed in latitude and season. 

The local dust storms which Viking detected in the southern hemisphere can be divided into 
two groups: those occurring within 10° to 15° ...of latitude of the receding south polar cap and those in 
the 10° to 30° S latitude zone. 

Common local dust storm localities in the southern hemisphere include: the region from Noctis 
Labyrinthus to Thaumasia Fossae, Argyre, Noachis-Hellespontus, Hellas, and the edge of the south 
polar cap. Common local dust storm sites in the northern hemisphere include: Chryse, Isidis-Syrtis 
Major, Cerberus, and Acidalia. Based on the number of observed clouds, an estimated 100 local 
storms occur annually. Most are associated with winds along the receding south polar cap and peak 
surface heating in the southern hemisphere near perihelion. Local dust storms form and dissipate in. a 
few days or less. A statistical treatment of the local storms is given in the "Occurrence Probability of 
Local Dust Storms” section. 

5. SURFACE ALBEDO 

The albedo of a surface is a measure of the incident solar radiation reflected by the surface. 

The albedo behavior of Mars was measured by the Viking IRTM instruments in the. range of 0.3 to 3.0 
pm [10]. Global dust storms on Mars play an important role in the deposition and removal of fine 
dust materials. The major variations in albedo were associated with the two global dust storms, of 
1977. During the seasons that are free of global dust storms the observed albedo was reasonably 
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uniform. Between the two storms, the atmosphere cleared substantially with many southern 
hemisphere regions returning to their prestorm albedos. The northern hemisphere atmosphere retained 
dust longer during the decay phase of the storms than the southern hemisphere. Southern hemisphere 
dark regions were not measurably brighter following the global storm suggesting little net deposition 
of dust (rate less than 0.15 to ~1.5 pm/year). In contrast, the northern hemisphere dark regions of 
Syrtis Major and Acidalia Planitia were measurably brighter following the storms, indicating a 
deposition of dust. These surfaces subsequently darkened over the following months returning to 
prestorm albedo values. Dust is also deposited in bright, low-inertia regions (at a rate of a few to ~45 
pm/year), where it remains. Remote sensing observations indicate that major dust deposits are 
currently located in three northern equatorial regions: Tharsis, Arabia and Elysium. They are covered 
by fine (~2 to 40 pm) bright (albedo > 0.27) particles. 

Considering solar collectors operation on the Martian surface in the context of dust deposition 
on the collectors and the change in the nearby surface albedo, we may take into account the following: 
the amount of solar radiation at perihelion and maximum wind velocities occur in the south causing 
storm initiation, whereas regional dust deposition occurs in the north suggesting net transport of dust 
from south to north. In the south, very small amounts of dust are deposited, whereas in the north 
there are regional dust deposits remaining on the surface for a length of time (few months), and there 
are regions where the dust accumulate and remain. 

The variation of the surface albedo on Mars with latitude and longitude is given in Table I for 
each 10° and is based on [11]. The albedo is taken into account also in solar radiation calculation on 
tilted surfaces. The mean albedo on Mars is about 0.27 and the maximum value is about 0.4. Table I 
may be used for clear days. For periods with dust storms we introduce an albedo function given by: 

al = maxja/j-, min(0.18t, 0.4)} (3) 

where al T is the albedo in Table I and 0.4 is the maximum albedo value. This function is based on 
the assumption that during dust storms the albedo increases above the value in Table I. It is also 
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assumed that dust storms that cause the atmospheric optical depth to increase above the value of 1.5 
will increase the surface albedo above its mean value of 0.27; thus obtaining the value of 0.18 
(0.27/1.5 = 0.18). A maximum value of 0.4 for the albedo is assumed regardless of the dust storm 
intensity. 

6. MARS SEASONS 

Mars has seasons comparable to those of Earth. However, the seasons are on the average about 
twice as long as on the Earth, corresponding to the greater length of the Martian year. Furthermore, 
they are distinctly unequal in duration as a result of the appreciable eccentricity of the Martian orbit. 
For that reason, the Martian year is not divided into months. Table II gives the duration of the 
Martian seasons in terrestrial and Martian days (a Martian day is defined as a sol, 1 sol = 24.65 hr). 
The Martian year is divided into degrees of areoeentric longitudes, corresponding to the position of 
Mars in its orbit around the Sun. Areoeentric longitudes L s = 0° and 180° correspond to the spring 
and fall equinox for the northern hemisphere, respectively, and L s = 90° and 270° correspond to 
northern and southern summer solstices, respectively. Figure 10 shows how many (terrestrial) days are 
in 1° of areoeentric longitude. This relationship is based on Mars orbit and the law of conservation of 
the angular momentum and is given by: 

, . 687 _q„ L f 2 )^ (4) 

360 [l + e cos (L s - 248°)] 

7. SOLAR RADIATION AT THE TOP OF MARS ATMOSPHERE 
The variation of the solar radiation at the top of the Mars atmosphere is governed by the 
location of Mars in its orbit and by the solar zenith angle. The solar radiation is direct beam radiation. 
The beam irradiance, in W/m 2 , is given by: 
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( 5 ) 



where L s is the areocentric longitude and 248° is the areocentric longitude of Mars perihelion. The 
Sun-Mars mean distance in astronomical units (AU) is 1.5236915; therefore, the mean beam irradiance 
at the top of the Martian atmosphere is: 1371/1.5236915 2 = 590 W/m 2 . The instantaneous beam 
irradiance is given by Eqs. (5) to (7). 



8. SOLAR RADIATION ON THE SURFACE OF MARS 
The variation of the solar radiation on the Martian surface is governed by three factors: (1) the 
Mars-Sun distance, (2) solar zenith angle, and (3) the opacity of the Martian atmosphere. The global 
solar irradiance G is the sum of the direct beam irradiance G b , the diffuse irradiance G d> and the 
albedo G^: 

' G = G b + Cr d + G al (9) 

The direct beam irradiance on the Martian surface normal to the solar rays is related by Beer’s 
law to the optical depth, t, of the intervening atmospheric haze: 
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G b = G ob ex P [-TW(z)l 


( 10 ) 


where m(z) is the air mass determined by the zenith angle z, and can be approximated, for zenith 

angles up to about 80°, by: 

m(z) * — — (11) 

cos z 

The zenith angle of the incident solar radiation is given by: 

cos z = sin <J> sin 5 + cos cj> cos 5 cos oo (1^) 

where 

<J> latitude 

5 declination angle 

co hour angle measured from the true noon westward 
The solar declination angle is given by: 

sin 5 = sin 5 0 sin L s (13) 

where 5 0 = 24.936° is the Mars Obliquity of rotation axis. The ratio of Mars to Earth length of day is 
24.65/24. It is convenient, for calculation purposes, to define a Mars hour, H, by dividing the Martian 
day (sol) into 24 hr. Using the same relationship between the Mars solar time T and the hour angle 
as for the Earth, we write: 

w = 15 T - 180 ( 14 ) 

This is shown in Fig. 11. The final solar radiation results can then be adjusted by the above 
ratio to correspond to actual (terrestrial) time. The sunset hour angle is given by [13]: 

co w = cos -1 (-tan <J) tan 5) (15) 

if |<t>| < tt/2 - 1 5 | . 
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If - tan 5 tan (J) > +1 the Sun will neither rise nor set for the day: 

or <|) < -Tt/l + 5 or <J> > jt/2 +5 polar night - > 

^ - If - tan 5 tan <j> = ±1 the Sun will be on the horizon for an instant 

only ( 16) 

If -tan 5 tan cj) < -1 the Sun will neither rise nor set for the day: 

or <|) > n/2 - 5 or <j) < n/1 - 5 polar day 

The number of Mars daylight hours is: 

T j - — cos -1 (-tan <f> tan 5) (17) 

d 15 , 

and the actual (terrestrial) number of daylight hours is obtained by multiplying Eq. (17) by the ratio 
24.65/24.' ' ' ' 1 ■ 

The solar irradiance on a horizontal Martian surface (G^ = 6) is given by: 

G h = G bh + G dh (l8) 

where 

G h global irradiance on a horizontal surface 
G bh direct beam irradiance on a horizontal surface 
G dh diffuse irradiance on a horizontal surface 

The diffuse irradiance of the Martian atmosphere may be a result of a mechanism different from 
that for the Earth atmosphere, nevertheless, to a first approximation, we will apply Eq. (18) as for 
Earth-terrestrial calculations. 

The global irradiance G h on a horizontal surface is given by: 

■ • ; • - v>t <h - e#«» t ( 19 ) 

where al is the albedo, and /( z/r,al) is called the normalized net flux function. The net solar flux 

integrated over the solar spectrum on the Martian surface was calculated by Pollack [11] based on 
multiple wavelength and multiple scattering of the solar radiation. Derived data from this calculation 
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are shown in Tables 111(a) and (b) for albedo = 0.1; and Tables IV(a) and (b) for albedo = 0.4. The 
parameters are the zenith angle z and the optical depth t. For albedos between 0.1 and 0.4 a linear 
interpolation can be used. Figures 12 and 13 describe the variation of the normalized net flux function 
with optical depth and Sun zenith angle for 0.1 and 0.4 albedo, respectively, for some data in 
Tables III and IV. The normalized net flux function can be approximated by a polynomial expression 
given by: 


5 5 i / \i 

f (z,T, al) = E £ E p(i,j,k) - V • | * (al) k 

_ /=0 j-0 k=0 100 J 

where p(i,j,k) are the coefficients of the polynomial given in Table V. These coefficients were 
obtained from fitting Eq. (20) to the data in Tables III and IV based on the criteria of least square 
error. The mean error is about 0.7 percent for the full range. For zenith angles up to 40° the error is 
much smaller. The largest error is for zenith angle of 80° and 85° and for T greater than 5. At these 
large angles and opacities, the error has a minor effect on the calculated daily insolation. The 
maximum error is about 7 percent. 

Using Eqs. (10) and (11), the beam irradiance G bh on a horizontal surface is given by: 



G bh - G ob cos z expj^Lj (21) 

and the diffuse irradiance is obtained by subtracting the beam from the global irradiances (Eqs. (19) 
and (21)). Figures 14 to 16 describe the variation of the global, beam and diffuse irradiances, 
respectively, on a horizontal Martian surface; and are given in pairs as functions of the optical depth T 
and zenith angle z. The irradiances were calculated based on Tables in and IV data, the mean 
irradiance of 590 W/m 2 and mean albedo of 0.27. The variation of the global irradiance G h , Eq. (19) 
is shown in Figs. 14(a) and (b). The beam irradiance G bh is obtained using Eq. (21) and is shown in 
Figs. 15(a) and (b). The beam irradiance shows a sharp decrease with increasing of the optical depth, 
and a relative moderate decrease with increasing of the zenith angle. The diffuse irradiance G dh is 
given in Figs. 16(a) and (b) and shows a sliding maximum with the variation of the zenith angle. 
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The solar radiation (global, beam and diffuse) variation (diurnally and daily) can be calculated 
based on the preceding equations and the /( z,T,al) function. For a given L s and (j), one can 
calculate the variation of the zenith angle z as function of the Mars solar time T using Eqs. (12) 
and (14). Referring to Figs. 3 and 4 for the given L g (or Eqs. (1) or (2)), the optical depth t is 
determined. With the albedo value Eq. (3), Tables III and IV (or Eq. (20)), and Eqs. (18), (19) 
and (21), one can calculate the solar radiation variation for the given day. The diurnal variation of the 
solar irradiances on a horizontal surface for L g = 141° and L g = 295° for Viking Lander location 
VL1 are shown in Figs. 17(a) and (b), respectively. Areocentric longitude L g = 141° corresponds to 
the lowest opacity of 0.35, and L g = 291° to the highest opacity of 3.6 (Fig. 3). The figures show 
clearly that at high opacity, the diffuse component dominates the solar radiation. The daily global 
insolation H h , beam isolation H b and diffuse isolation H d can be calculated by integrating the 
iiradiances over the period from sunrise to sunset. Figures 18(a) and (b) shows the daily global 
insolation H h on a horizontal surface, in kWhr/m 2 -day at VL1 and VL2 over a Martian year. The 
albedo at VL1 is 0.22 and at VL2 is 0.24. Also shown is the percentage of the diffuse insolation 

on a horizontal surface. During the spring and summer where the optical depth is relatively 
low, the diffuse insolation comprises about 50 percent of the global insolation, whereas for autumn 
and winter, the optical depth is high and therefore the diffuse insolation is about 90 percent. 

A different view on the amount of solar radiation on Mars is the cumulative irradiance for the 
entire year. This is presented in Fig. 19. The ordinate is the instantaneous irradiance and the abscissa 
is the percent of yearly sunshine hours. The figure shows, for example, that during about 55 percent 
of the yearly sunshine hours the irradiance on a horizontal surface at VLJ is above 200 W/m . 

The variation of the daily average irradiance in W/m and the daily insolation in Whr/m -day 
ottsa horizontal surface for the entire planet (for latitudes other than for the Viking Landers) will now 
be given based on the models of the optical depth functions, t(<J),L s ), Eqs. (1) and (2). The average 
irradiance is obtained by dividing the insolation value by the number of daylight hours for the given 
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day, Eq. (17). The variation of the global irradiance, G h , with latitude <|> and areocentric longitude 
L $ on a horizontal surface using Eq. (1) is shown in Fig. 20, and the daily global insolation, H^, is 
shown in Fig. 21. The variation of G h and H h on a horizontal surface using Eq. (2) is shown in 
Figs. 22 and 23, respectively. 

The two global dust storms that were observed in 1977 may be considered as one of the worst 
years of dust activities on Mars. However, there are many Mars years in which no global storm 
occurred. In Fig. 24 we compare the daily global insolation at VL1 and VL2 for the observed 
opacities in 1977 (see Figs. 3 and 4) and for an assumed Martian year with no global dust storms and 
with a constant optical depth of T = 0.5. 

9. OCCURRENCE PROBABILITY OF LOCAL DUST STORMS 

Local dust storms may be related to clouds and hazes of the Mars atmosphere. The occurrence 
of the local dust storms is of statistical nature and, therefore, has to be treated accordingly. The 
estimation of the occurrence probability of these storms at a given location on Mars is based on 
Table 3 of [14]. The clouds and hazes in percentages are given in [14] for 6 latitude and 10 seasonal 
bins. The hazes are divided into three categories according to the visibility of surface features through 
the obscuration. "Thin” haze (visible optical depth T < 1); "moderate” haze (t » 1); and "thick” 
haze (T > 1). Each hemisphere is divided into three latitude bins: equatorial, 0° to 30° latitude, 
mid-latitude, 30° to 60°; and polar, 60° to 90°. The season bins, include 0° to 45°, 45° to 80°, 

80° to 125°, 125° to 160°, 160° to 200°, 200° to 245°, 245° to 270°, 270° to 290°, 290° to 325°, and 
325° to 360° of areocentric longitude L s . The occurrence probability P was defined as the number 
of clouds or hazes divided by the number of observations for a given period of time. Based on this 
information we developed a statistical model for the occurrence of the local dust storms. All three 
types of hazes were combined into a single probability number. The clouds were not considered since 
they are usually smaller in spatial extent. The occurrence probability in [14] reports on clouds and 
hazes seen by Viking observers in a given bin area.. Our purpose is to translate these events into an 
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occurrence probability as seen by an observer on Mars surface. Since the area. of a bin is large, a 
cloud occurrence in this bin may not be detected by an observer on the surface. In order to estimate 
the occurrence probability of a local dust storm affecting an area on Mars surface, the probability P 
has to be normalized with respect to areas, i.e.. 



where A is the area covered by a local dust storm, and A„ is the bin area given by: 

A m = 27ri? 2 (sin4> 1 - sin<j> 2 ) ( 23 ) 

where R is the Mars radius, and 4> 1 and <j) 2 are the latitudes defining the bin. Table VI gives the 
occurrence probability P r (in percentage) for a local dust storm covering an area A of 10 5 km 2 for 
the given latitudinal and seasonal bins. 

9.1 Occurrence Probability of a Number of Local Dust Storms 

It is possible to estimate the occurrence probability P(N) pf local dust stontts for a given length 
of time t at a given location. This is given by: 

e[/>(W)] = P r L (24) 

f s 

where E[P(N)] is the expected value of the occurrence for local dust storms and tg is their duration. 
To calculate the occurrence of N storms, P(N), we assume a Poisson distribution where an 
occurrence of a storm is independent on the occurrence of subsequent storms, i.e., 


P(N)=k N (25) 

For an observation time t equals to the duration time t g of the local storm it is possible to show that 


e[P(N)} = E(k N ) = P r = 


( 1 \2 




(26) 


or 


k = 


_ 2 P r + 1 - ftp,. + 1 


(27) 


2Pr 
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where k is the occurrence probability for one local storm in the period t = t s . Small values of P r 
(Eq. (27)) yields k * P r , i.e., the occurrence probability of one local dust storm is given by P r in 
Table VI. The greatest probability as indicated in (Table 3, [14]) is for the latitudinal bin -90° 
to -60° and seasonal bin 245° to 270°. According to Table VI, P r = 0.06852 and from Eq. (24) 
k = 0.0607, i.e., k » P f . The occurrence probability of a local dust storm is (Eq. (25)) P( l) = k, i.e., 
6.07 percent. The occurrence probability of two storms is P(2) = k 2 = 0.0037. 

Local dust storms may be taken into account in solar radiation calculation. The occurrence 
probability P r represents the local dust storms in the solar radiation calculation. It is used for the 
optical depth models for latitudes other than for the measured values at the Viking Landers VL1 and 
VL2. This is because the measured optical depths already contain the effect of the local storms, if 
they really occurred. To determine whether a local dust storm occurred for a given latitude <]) and 
day L s , a number is randomly selected in the range between 0.0 and 1.0. If the value P r in 
Table VI is greater than the randomly selected number we assume that there is a local dust storm. The 
models for the optical depths assumes t = 0.5 for clear skies. The optical depth for a local storm is 
taken as t = 1 . The effect of the local dust storms on the daily global insolation 11^ on a horizontal 
surface is calculated according to: 

H h = H h\ 0 ~ P r) + H h\ p r (28) 

T-0.5 T-1.0 

where P r is given in Table VI. 

An example of this calculation is shown in Fig. 25. The latitude cj> = -70° (i.e., very far south) 
was selected for the calculation because this region is most affected by local storms. In order to 
exaggerate the effect, the dust storm probability was arbitrarily doubled over that in Table VI. The 
assumed optical depth t was 0.5 and for a year without global dust storms. For comparison purpose 
Fig. 25 also shows the curve for t = 1.0, showing the insolation during a local dust storm. The effect 
of the local dust storms is noticed by comparing % for T = 0.5 with the H h . The largest effect, in 
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this example, is 1.3 percent for L s “ 265. On the average, local dust storms have a minor effect on 
the daily global insolation. 

9.2 Number of Days with a Local Dust Storm in a Year at a Given Location 

Local dust storms may be taken into account in sizing the capacity of the energy storage. Days 
with local dust storms may be considered as days with low irradiance or "dark days.” The number of 
days, in units of L s , in a Martian year with local dust storms may be determined for each latitudinal 
bin A(j> by: 

d ts = E P,.(A<|), A L s ) x A L s (29) 

AL, 

where AL S is the seasonal bin. This is presented in Table VII. 

If the duration of a local dust storm is 5° of L g , the number of local dust storms for each 
latitudinal bin per year is obtained by dividing the values in Table VII by 5°. 

CONCLUSIONS 

The solar radiation model for Mars has been refined to include an extension of the net solar flux 
function /(z,r,al), surface albedo, optical depth function for different latitudes and a statistical 
treatment of local dust storms. The paper presents a procedure and solar radiation related data from 
which the diurnally and daily variation of solar radiation on Mars were calculated. 

The two global dust storms that were observed in 1977 make it one of the worst years of dust 
activities. Even so, the opacities are relatively low during the northern spring and summer for more 
than a half Martian year (or more than a terrestrial year) for which the insolation is relatively high 
varying from about 2 to 4 kWhr/m 2 -day. These seasons are the likely ones during which manned 
missions may be carried out. It should also be noted that there are many Mars years in which no 
planet-encircling dust storm has occurred; less frequently, there may be years in which not even 
regional storms occur. One of the most important results of this study is that there is a large diffuse 
component of the solar insolation, even at high opacity, so that solar energy system operation is still 
possible. 
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The insolation model for Mars presented in this paper enables more accurate engineering studies 
and estimates of solar cell array size and performance than has heretofore been possible. 
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NOMENCLATURE 

A area covered by a local storm 

al Mars surface albedo 

d number of terrestrial days 

d & number of days (in units of L s ) of local dust storms 

G b beam irradiance on Mars surface 

G h ,G b h»G d h global, beam and diffuse irradiance on Mars horizontal surface 
G ob beam irradiance at the top of Mars atmosphere 

H h ,H b h>Hdh global, beam and diffuse daily insolation on Mars horizontal surface 
Subscripts: 

b * direct (beam) values 

d diffuse values . 

h horizontal values 

o values on top of Mars atmosphere 

Other values: 

e eccentricity = 0.093377 

k occurrence probability for one local storm 
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L s areocentric longitude (position of Mars in orbit around the Sun) 

m(z) air mass 

r Sun-Mars distance 

R Mars radius ' "" 

S solar constant =1371 W/m 2 at the mean Sun-Earth distance of 1 astronomical unit (AU) 

T Mars solar time 

t time 

t s duration of a local storm 

P r occurrence probability of a local storm 

P(N) occurrence probability of local storms 

T d Mars daylight hours 

z zenith angle 

5 declination angle 

5 0 obliquity 

0 true anomaly 

r optical depth - 

<j> latitude 

co hour angle 

co ss sunset hour angle 
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TABLE II. - MARS SEASONS 


Areocentric 

longitude 

Season 

Duration of the season 
Mars 

or me oun, 

L s> 

deg 

Northern 

hemisphere 

Southern 

hemisphere 

Martian 

days 

Terrestrial 

days 

0 to 90 

Spring 

Autumn 

WM 

. 199 

90 to 180 

Summer 

Winter 

wSm 

182 

180 to 270 

Autumn 

Spring 

mm 

146 

270 to 360 

Winter 

Summer 

■gs 

160 

687 
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TABLE 111(a). - NORMALIZED NET FLUX FUNCTION /( z,T,al) 
AT THE MARTIAN SURFACE 





















TABLE 111(b). - NORMALIZED NET FLUX FUNCTION /(z,T,al) 
AT THE MARTIAN SURFACE 
Albedo 0.1 
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TABLE IV(a). * NORMALIZED NET FLUX FUNCTION /(z,T,al) 
AT THE MARTIAN SURFACE 


Albedo 0.4 


Zenith angle z, deg 


iepth, 

r 

0 

5 

10 

15 

0.10 

0.596 

0.595 

0.594 

0.593 

.15 

.592 

.590 

.590 

.589 

.20 

.587 

.586 

.585 

.584 

>25 

.583 

.582 

.581 

.579 

.30 

.578 

.577 

.576 

.575 

.35 

.574 

.573 

.572 

.570 

.40 

.569 

.568 

.567 

.565 

.45 

.565 

.564 

.563 

.560 

.50 

.560 

.559 

.558 

.556 

.55 

.555 

.555 

.553 

.551 

,60 

.550 

.549 

.548 

.546 

.65 

.545 

,544 

.543 

.541 

.70 

.540 

.540 

.538 

.536 

.75 

.535 

.535 

.533 

.531 

.80 

.531 

.530 

.528 

.526 

.85 

.526 

.525 

.523 

.520 

.90 

.521 

.520 

.518 

,515 

.95 

.516 

.515 

.514 

.510 

1.00 

.511 

.511 

.509 

.506 

1.05 

.506 

.506 

.504 

.500 

1.10 

.501 

.501 

.499 

.495 

1.15 

.497 

.496 

.494 

.490 

1.20 

.492 

.491 

.489 

,485 

1.25 

.487 

.486 

.484 

.480 

1.30 

.482 

.481 

.479 

.475 

1.35 

All 

.476 

.474 

.470 

1.40 

All 

.471 

.469 

.465 

1.45 

.467 

.467 

.465 

.460 

1.50 

.463 

.462 

.460 

.456 

1.55 

.458 

.457 

.455 

.451 

1.60 

.453 

.452 

.450 

.446 

1.65 

.449 

.448 

.445 

.441 

1.70 

.444 | 

,443 

.440 

.436 

1.75 

.440 ! 

.438 

.435 

.431 

1.80 

.435 

.434 

.431 

.427 

1.85 

.431 

.429 

.427 

.422 

1.90 

.426 

.425 

.422 

.417 

1.95 

.422 

.420 

.418 

.413 

2.00 

.417 

.416 

.413 

.408 

2.10 

.408 

.407 

.404 

.399 

2.20 

.399 

.398 

.395 

.391 

2.30 

.391 

.390 

.387 

.382 

2.40 

.383 

.382 

.379 

.374 

2.50 

,375 

.374 

.371 

.367 

2.60 

,367 

.365 

.363 

.358 

2.70 

.358 

.357 

. .354 

.350 

2.80 

.350 

.349 

.346 

.342 

2.90 

.343 

.342 

.339 

.334 

3.00 

.336 

.335 i 

.332 

.327 

3.20 

.322 

.321 

.318 

.314 

3 . 40 , 

.308 

.307 

.304 

.300 

3.60 

.294 

.293 

.290 

.286 

3.80 

.282 

.281 

.278 

.274 

4.00 

.271 

.270 

261 

.263 

4.50 

.244 

.243 

.241 

.237 

5.00 

.219 

.218 

.216 

.212 

5.50 

.196 

.196 

>194 

*190 

6.00 

.176 

.175 

.173 

.170 


20 

25 

30 

35 

40 

0.592 

0.591 

0.590 

0.589 

: 0.588 

.588 

.586 

.585 

.583 

.582 

.583 

.581 

.579 

.577 

.575 

.578 

.576 

.574 

.571 

.569 

.573 

.571 

.568 

.565 

.562 

.568 

.565 

.563 

.559 

.556 

.563 

.560 

.557 

..553 

.549 

.558 

.555 

.552 

.547 

.542 

.553 

.550 

.546 

.541 

.535 

.548 

.544 

.540 

.535 

.529 

.543 

,539 

.534 

.529 

.522 

.538 

.533 

.528 

.522 

.516 

.532 

.527 

.522 

.516 

.509 

.527 

.522 

.516 

.510 

.503 

.522 

.517 

.511 

.504 

.496 


l J 


















TABLE IV(b). - NORMALIZED NET FLUX FUNCTION /( z,T,al) 

AT THE MARTIAN SURFACE 

Albedo 0.4 
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TABLE V. - NORMALIZED NET FLUX FUNCTION COEFFICIENTS, /(z,r,al) 


k = 0 

i 

j \ 

0 

1 

2 

3 

4 

5 

0 

1.002800 

-0.228681 

0.019613 

0.000231 

-0.000130 

0.000003 

1 

-.450073 

1.335955 

1 .131691 

.402126 

-.063967 

.003758 

2 

5.566705 

-16.912405 

13.739701 

-4.756079 

.743740 

-.043159 

3 

-22.471579 

64.909973 

-52.509470 

17.997548 

-2.786548 

.160340 

4 

36.334497 

-101.800319 

79.895539 

-26.762885 

4.0741 17 

-.231476 

5 

-20.420490 

53.207148 

-39.949537 

12.977108 

-1.931 169 

.107837 

k = 1 ; 

■ 

0 

1 

2 

3 

4 

5 

0 

0.009814 

0.226139 

-0.1 17733 

0.030579 

-0.004090 

0.000218 

1 

-.156701 

.396821 

-.313648 

.099227 

-.013508 

.000651 

2 

1.361 122 

-3.758111 

3.007907 

-.987457 

.141693 

• . -.007320 

3 

-4.365924 

12.539251 

-10.394165 

3:486452 

-.513123 

.027401 

4 

5.991693 

-17.498138 

14.291370 

-4.765323 

.703675 

-.037960 

5 

-2.915099 

8.275686 

-6.59312 5 

2.173999 

-.320308 

.017335 


TABLE VI. - PERCENT OCCURRENCE PROBABILITY OF A LOCAL DUST STORM 


Latitude, 

deg 

L s , deg 

0 to 45 

45 to 80 



160 to 200 

200 to 245 

245 to 270 

270 to 290 

290 to 325 

325 to 360 

-90 to -60 

2.816 

i.583 


0.000 

3.079 

6.166 

6.852 

4.111 

4.756 

1.541 

-60 to -30 

.706 

.399 

.197 

.928 

1.368 

1.968 

1.754 

1.289 

1.426 

.418 

-30 to 0 

.262 

.703 

.323 

.227 

.359 

.594 

1.156 

.648 

.966 

.383 

0 to 30 

.238 

.703 

.645 

.620 

.497 

.514 

.774 

.433 

.716 

.414 

30 to 60 

.527 

.654 

.806 

1.022 

1.044 

1.177 

1.304 

1.756 

.967 

.903 

60 to 90 

3.494 

2.401 . 

2.372 

2.668 

4.275 

5.483 

.000 

.000 

.000 

.000 


Table vii. - number of 

DAYS WITH A LOCAL 
DUST STORM* IN A 


MARTIAN YEAR 


Latitude bin, 
deg 

Days, L s , 
deg 

-90 to -60 

10.57 

-60 to -30 

3.65 

-30 to 0 

1.89 

0 to 30 

1.97 

30 to 60 

3.47 

60 to 90 

8.59 
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Optica! depth, t Optical depth, t 



Areocentric longitude, L s , deg 

Fig. 1 - Measured and discretized optical depth at Viking Lander VL1 and years of great dust storms. 
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Areocentric longitude, L s , deg 

Fig. 2 - Measured and discretized optical depth at Viking Lander VL2 
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OPTICAL DEPTH 




Fig. 3 - Optical depth at Viking Lander VL1 - in detail. 


31 


OPTICAL DEPTH, 
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AREOCENTRIC LONGITUDE, L s , DEG. 

Fig. 4 - Optical depth at Viking Lander VL2 - in detail. 
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0 30 60 90 120 ISO ISO 210 240 270 300 330 360 

AREOCENTRIC LONGITUDE, L s , DEG. 

Fig. 5- Variation of optical depth with areocentric longitude for latitudes of the Viking 
Landers VL1 (22.3°N) and VL2 (47.7°N), and for the latitude -BO® of the two global 
storm sources. * : ' 
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Areocentric Longitude. LS. Deg. 

Fig. 10 - Earth days per one degree of areocentric longitude. 



T = 12:00 


Fig. 1 1 - Solar time and hour angle relation. 
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Surface albedo = 0.1 
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(a) Effect of optical depth with Sun zenith angle as a parameter 
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(b) Effect of Sun zenith angle with optical depth as a parameter 

Fig. 16 - Variation of diffuse Irradiance with optical depth and sun zenith angle, 
on a horizontal surface, at mean irradiance of 590 W/m 2 
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AREOCENTRIC LONGITUDE, Ls , DEG. 


Fig. 23 - Variation of daily global insolation H h (Whr/m 2 day) with latitude and arcoccntric longitude on a horizontal 
surface of Mars based on the optical depth function cq. (2). 




0 90 180 270 360 


Areocentric longitude, L s , deg 

Fig. 24 - Daily global insolation at VL1 and VL2 for the 
observed opacities at VL1 and VL2 and for t = 0.5. 
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Fig. 25 - The effect of local dust storms on the daily global insolation on a horizontal surface 
at latitude of 70° S. 
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